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This thesis includes a survey of the literature through 
1967 on the methods of determination of the thermodynamic proper- 
ties of nitrogen-oxygen mixtures in the tvjo-phase, liquid-vapor 
region. The virial expansion, empirical equations of state, pseu- 
docritical constants, excess functions, latent heats, and graphi- 
cal methods are considered. 

The methods just enumerated are evaluated for their simpli- 
city, feasibility, and accuracy. On the basis of the state of the 
art and the* availability of appropriate data, thermodynamic dia- 
grams are constructed for 1 , 5 i 10 and 20 atmospheres of pressure 
in the tvro-phase region. Diagram coordinates are enthalpy- 
composition and entropy-composition. The diagrams at atmospheric 
pressure show good correlation with knovm data, ‘while those at 
higher pressures are best estimates using excess functions and 
idealizing assumptions. 

Thesis Supervisor: J.L. Smith, Jr. 

Title: Associate Professor of Mechanical Engineering 






, 






;'$l * 



4;|4 



Title page 
Abstract 



TABLE OF 'CONTENTS 



Table of Contents 
List of Symbols 

I, Introduction 
II, Procedure 
III, Results 



1 , 


Phase Equilibria 


2 , 


The Virial Equation 


3 . 


Pseudocritical Properties 


4 , 


Empirical Equations of State 


5 . 


Partial Properties 


6 , 


Excess Functions 


7 . 


Calorimetry 


8 , 


Latent Heat 


9 . 


Graphical Methods 


10 , 


The Thermodynamic Properties 



IV, Discussion of Results 
V, Conclusions 
VI, Recommendations 
VII, Appendices 



A, 


Equilibrium Diagrams 


B, 


Numerical Technique of Barker 


C, 


Selected Data 



111 



LIST OF SYMBOLS 



a, b; a*, b^; a'^ , bY single substance, mixtvire and pure com- 

ponent of mixture constants in Redlich- 
Kv:ong eqn, of state, 

a, A, b, B, c, C constants in Barker's expression for 

B(T), C(T) 2nd & 3rd virial coefficients for a single sub- 

stance , 

Bj^., B(T,x), Bjjj 

2nd & 3rd virial coefficients for a mixture, 

jk‘ C(T,x) , 

®AA* ®bB’ ^AAA* ^bBB ^ virial coefficients for pure 

components in a binary mixture. 



®AB» ^ABB* *"AAB ^ interaction or cross virial coeffi 

cients for a binary mixture. 



«• * 

B(T), C(T) 2nd & 3rd universal, single substance, virial 

coefficients , 

"hard sphere" molecular volume =(2/3)7T 

C® specific heat at constant pressure for perfect gas or at 

^ zero pressure, 

D(T,x) etc, — •- 4th virial coefficient, 

E internal energy, 

f fugacity of a single substance, 

fugacity of component i in solution 

fj^ fugacity of pure component i ?.t same T,P of the mixture, 

G Gibbs free energy = H - TS, 



H enthalpy = S + PV 

H, » /\h , differential heat of vaporization and condensa- 

^ tion at constant temperature. 



absolute joules 

k Boltzmann constant = 1,38054 x 10^^ 

L, M, N — defined by Eq. (58). 
n number -of moles. 



iv 
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P pressure. 

pure-component pressure of component i In a gas mixture, 

partial pressure of component i in ideal gas mixture, 

= x|^ P (Dalton’s Lavj) , 

R universal gas constant = 0.0820535 liter atm/mole °K 

R any extensive property 

partial property of component i in solution, 

S entropy 

T absolute temperature, 

base or reference temperature. 

U intermolecular interaction energy, 

V --- molal or specific volume 

partial molal volume of component i in solution, 

mole fraction of a component i in a mixture = n^/n, 

Z compressibility factor = PV/RT 

superscripts 
E excess property, 

1 property of the liquid phase; property of the liquid line, 

V property of the vapor phase; property of the vapor line. 

m property of mixing, 

0 zero pressure state or perfect gas state at any pressure. 

subscripts 

b refers to the normal boiling point. 

1 pure component property of a mixture. 

m property of a mixture. 

c critical property, 

cm pseudocritical property. 
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r reduced state. 

s reference state. 

Greek symbols 

CX volatility ratio. 

activity coefficient. 

Q general function of temperature. 

£ characteristic molecular energy. 

(j characteristic molecular distance, 

finite change of quantity, 

yU chemical potential = G 

All state properties are given in their molal or specific proper 
ty form unless otherv:lso noted. 
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INTRODUCTION 



It is the pui'pose of this thesis to present the results of a 
search of the literature and the subsequent analysis of the find- 
ings of that search. The subject of the inquiry was the thermo- 
dynamic properties of nitrogen-oxygen mixtures in the two-phase, 
liquid-vapor region. The reasons for the study of the elemental 
properties of any basic substances are self-evident especially 
for substances as fundamentally important as the mixtures of nitro- 
gen and oxygen. It should suffice to say that the hunt carried 
the author through the fields of cryogenics, low temperature phy- 
sics, high pressure physics, thermodynamics, physical chemistry, 
molecular chemistry and statistical mechanics. It is the intent 
of this thesis to extract from these areas of science and engi- 
neering the most appropriate analytical, graphical and machine 
methods of determination now available and to present the proper- 
ties determined therefrom, 

Acknov^ledgement and appreciation are expressed for the 
aid, assistance, and direction of Professor J,L, Smith, Jr, 
of the MIT Cryogenics Laboratory, 
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PROCEDURE 

The direction of this work consisted of two distinct phases. 
The first part of the effort was directed toward a thorough search 
of the scientific and engineering literature. This was accom- 
plished by a personal search of the Chemical Abstracts through 
1967. Machine programmed searches v^ere conducted for the author 
by the Cryogenics Engineering Laboratory of the National Bureau of 
Standards and the Defense Documentation Center. These initial 
efforts were principally centered on nitrogen-oxygen mixtures and 
the thermodynamic approaches to the derivation of properties of 
binary mixtures. As the work progressed, molecular theory became 
of interest, and the extensive bibliographies in the major texts 
of the field vrere utilized. At the end of each subsection of the 
Results, the references are listed vrhich deal explicitly vrith 
that subsection as annotated in the text. 

The second part of the program vras the analysis of the data 
which had been found. It must be said that the determination of 
the properties of a mixture would be a relatively easy problem if 
the mixture behaved simply as the average of its constituents. 
Unfortunately, not only are the constituents themselves not ideal, 
but also they interact vrith each other in a non-ideal manner. 

It is intended in this thesis to discuss the PVT behavior of 
mixtures 8 .nd the various methods of deriving their thermodynamic 
properties: enthalpy and entropy. Specifically the next section 

will cover the follovring: 

a, the theory of corresponding states and pseudocritical con- 



stants 
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b, empirical equations of state, 

c, the virial expansion, 

d, partial property correlations and mole fraction averaging, 

e, excess functions as deviations from ideality, 

f, graphical derivations, 

g, latent heat techniques. 

No attempt has been made to make interpretations based on statis- 
tical mechanics and its application to molecular theory. The read- 
er interested in approaches based on this field must consult the 
vjorks of Eyring, Kirkwood, Lewinson, Hirschf elder et al for their 
treatment on the formal and model approaches in this area. 

The final section of the Results vrill correlate the discussed 
data and provide a graphical presentation of the enthalpy and en- 
tropy of binary mixtures of nitrogen and oxygen at various pres- 
sures , 
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RESULTS 



Phase Eouilibria 

Perhaps the information most consistently required in calcu- 
lations concerning the tv;o-phase region is the equilibrium behav- 
ior of the system as defined by Gibbs [l]« The study of vapor- 
liquid equilibrium would obviously be simplified if the phases 
behaved as ideal solutions, for then vre should only require infor- 
mation about the pure components in order to develop values for 
the mixtures , 

Deviations from ideal behavior occur in both the liquid and 
the vapor phases as a result of the forces that exist between the 
molecules. The magnitude of the intermolecular forces depends on 
the nature of the molecules and the average distance betv;een them 
as expressed by 



This expression of the principle of corresponding states 

says that the interaction energy of a pair of molecules A and B 

at a distance r apart may be expressed in terms of a universal 



characteristic energy and distance respectively for the pair of 
molecules. In the liquid state the molecules are relatively 
close together, and the force fields are strong, so that differ- 
ences in shape and charge distribution vrlll cause deviations in 
behavior. In the gas phase the force fields are not as strong, 
and at lov; pressures ideal mixture behavior is a good approxima- 
tion, Mote that a perfect gas is defined in terms of negligible 






( 1 ) 




function f for all A and B, The parameters £ and are a 



I 




5 



forces between molecules, vj’hereas an ideal mixture or solution 
requires only the same forces between like and unlike molecules. 
Hence, at lov; pressures ideal solution deviations are principally 
due to the liquid phase vihile at high pressures both vapor and 
liquid phases contribute. 

It is perhaps instructive to review the phase behavior of 
mixtures under pressure [4,5,6], In Fig. 1 the curves OCj^g 

represent the usual vapor pressure curves for the pure com- 
ponents terminating at their respective critical points. For the 
binary system, however, the phase boundary is the region enclosed 




Figure 1. Phase boundaries of pure components and mixtures. 
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by the curve within which both the liquid and the vapor coex- 

ist. The critical point of the fixture, C , is defined as the 

ra 

state at which the properties of the vapor and liquid phases in 
equilibrium are Identical, If the phase boundary curves of sever- 
al mixtures of different compositions are drawn, the locus of 
their critical points form an envelope indicated by the dashed 
line in Fig. 1. The maximum pressure and temperature for the ex- 
istence of two phases are called the cr icondenbar , Cp, and the cri- 
condentherm, Cp , respectively [?j. These are illustrated in Fig, 

2 as is the phenomena of retrograde condensation 8 , 





Figure 2 



Figure 3 



In Fig, 2 it can be seen that if a line of constant temperature 
is follov7od, condensation v/ill begin at point a and will continue 
v7ith decreasing pressure until maximum condensation is reached at 
point b, at vrhich tine vaporization ’will occur until the system is 
again all vapor at point c. A similar phemomcna vrill take place 
at constant pressure with vaporization occurring with decreasing 



I- 
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temperature from points d to e, • 

The diagrams raost frequently used in representing equilibri- 
um data are those of pressure or temperature as a function of 
molar composition. Fig. 3 shov:s the trend of the phase boundary 
curve V7ith Increasing pressure when the mixture critical tempera- 
ture and pressure fall betvxecn those of the pure components. The 
area of two-phase coexistence decreases with increasing pressure. 
The phase boundary curve breaks avray from the puj:e component vrhose 
critical temperature is lox^rest and recedes along the critical en- 
velope, finally disappearing altogether at the higher critical 
point. 

Turning to nitrogen-oxygen mixtures in pa.rticular, the earli- 
est investigation of their equilibria was conducted by Baly in 
1900 at atmospheric pressure |^9j . The temperature scale used by 
him v:as in error, and his vapor composition measurements proved 
to be innaccurate, Inglis [lo] la I 906 mads measurements at 
74,7 and 79.07°K resulting in vapor pressures betxveen 1 and 32 
atmospheres. From this work an empirical formulation was made 
for th$ tabulation of isotherms and isobars. 

During the next few decades the only measurements vjere those 
of Trapeznikova and Shubnikov [ 13 ] ia 193^+ oa an isotherm of 
85°K, and the static measurements of Sagenkahn and Fink [l^] 

1 and 2 atm in 19^^. 

Little else was done until 1955 xvhen Armstrong and his col- 
leagues at the National Bureau of Standards used a circulation 
method to make accurate determinations of the compositions and 
vapor pressures at 65 1 70, and 77.5°K [l5j» Tvto years later 
Cockett conducted extensive measurements at 1.3158 atmospheres 



^l6j ’.'.'ith 'results reasonably close to those of Dodp:e and Dunbar. 

Recent significant measurements are those of Din who 

used the static method for nine intermediate mixtures at pressures 
to 10 atm. His results likevrise shovr those of Dodge and Dunbar 
to be reasonable as do the experiments of Latimer in 1957. 

In 1961 Yendall and Olzewski follovjed with a report on the satur- 
ation properties of nitrogen-oxygen mixtures [l9j . 

The V/ADD Technical Report of 19^1 [^oj presents a compila- 
tion of the data of Armstrong et al, Dodge and Dun"bar and Sagen- 
kahn and Fink on T-x diagrams. The data of Sagenkahn and Fink at 
1 and 2 atm and of Cockett at 1.3158 atm are presented on P-x 
diagrams. The graphics use a base of Hilsenrath’s [ 21 J pure com- 
ponent vapor pressure values. Tables of selected experimental 
values of all of the above are given. 

Air Product and Chemicals, Inc. also published under govern- 
ment contract a technical report [rrJ with a P-T diagram of vapor 
and liquid c.urves based on the data of Dodge and Dunbar for the 
tv;o-phase region and the vapor pressure data of Hoge for oxygen 
1 ^ 23 ^ and -Michels et al, Friedman and VJhite, Dodge and Davis, and 
Henning for nitrogen ^24,25,26,27]. The critical envelope is 
based on the vrork of Kuenen et al j^28,29] and can be considered 
precise to 't2°K, The vrork of Kuenen et al also Indicates the 

I 

presence of retrograde condensation over small ranges near the 



critical points. 

The most recent vrork is that of Haselden and Thorogood in 
1963 at high oxygen concentrations [ 30 ] which shovred a good corre- 
lation with the earlier vrorks of Din and [^atim.er belovr 96 ^mole 
oxygen. They demonstrated, hovrever, that the relative volatility 
falls off at high oxygen concentrations. This is alleged to be 
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indicative of the occurrence of association, that is, the vreak 
chemical interaction between the oxygen and nitrogen molecules. 

For the vrork at hand, reliance has been placed primarily in 
the data of Dodge and Dunbar, because it is considered to be con- 
sistent and detailed over a wide range of values. Their data 
can be considered accurate to l0,l°K at pressures near atmospheric 
and to about 1]1°K at pressures above seven atmospheres [szj , 
Saturation values for the pure components have been taken for 
this vrork from Stevrart for oxygen and Strobridge for nitro- 



gen |^ 32 j , since the data of these investigators vras used for all 
other pure component property values, Stewart’s empirical vapor 
pressure equation v;as fitted to the data of Hoge [ 23 ] vrhile that 
of Strobridge fits the data of Friedman and White and Arm- 



Individual temperature-composition diagrams have been drawn 
for 1, 5 1 10 20 atm and are presented in Appendix A, 
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The Yirirfl Equation 



The virial equation of state is simply a series expansion of 
the compressibility factor in terms of density. 

Z = Py = 1 4 . B(TjX) ^ C(T,x) , D(T,x) ^ (2) 

RT V ^ ■ v3 ■ ^ * 

Because of the molecular interactions, the virial coefficients are 
functions of both temperature and composition and can be vjritten 
as [l] 

n n 

B = B(T,x) =E E ^i^ 1 ®ii and ( 3 ) 

“ i=.-A j=A ^ ^ 

n n n 

(It) 

^i=A j=A k=A 

For a binary mixture these become 

= "^A^aa + .2XaXbBab + 

“ ^ppAAA ^^A^B^AAB ^^a’^B^BBA ^B^BBB 

vxhere x^ is the mole fraction of component A, 

Molecular theory tells us that for non-polar molecules (angle 
independent molecules possessing no permanent charge distribution) 
such as those being investigated, calculation of their virial co- 
efficients depends on the cho 3 ce of proper potential functions. 

The potential function represents the attractive and repulsive 
forces betvreen molecules and many are found in the literature 
[ 2 , 3 I . Although these forces are generally exponential in ch^r- 



acter, they have generally been simplified to reciprocal povjcr 
forms . 

The most popular relationship is of the form 
U(r) = K^r“^- K2r‘® (7) 

where the first term represents the repulsive potential energy . 
and the second that of the attractive contribution. Where n and 
m are chosen as 12 and 6 respectively there results the well knovm 
Lennard-Jones potential (compare vrith Eq, 1 ), 

U(r) = {(j/r)^^^{(j/r)^ (8) 

It can be shown that by using the Lennard-Jones potential the 

pure component virial coefficients can be expressed as 

B(T) = f(T",b ) or = B(T)/b = fCT^") ( 9 ) 

o o 

vrhere b^ is the "hard sphere" volume defined by (2/3)/T-”^^ ^nd 
= T/(C/k). Similarly 

c''(T) = C(T)/bQ^ = f(T''/6'/k) • ( 10 ) 

Eqs, ( 9 ) & ( 10 ) have both been tabulated and plotted. 

The interaction or cross virial coefficients are signifi- 
cantly more difficult to determine. For an ideal solution obey- 
ing Ar.agat’s Rule 

S.AB^^ ^ (11) 
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v.’hich is the Levjis-Randall rule for mixtures of simple gases and 
results in the following for a binary mixture: 



(12) 



Cm=^A^A AA'’ ^ %A ~^BB ^ 

It has been shox^n » however, that as 
universal function of kT/^^^, so too is ^ab^"^ ’^^'^(j'ab^ ^ function 
of kT/ ^ The measurement of second virials in mixtures of di- 

lute gases where the third virials are negligible can be used to 
verify the various combining rules derived on the basis of this 
relation. The most popular and sim.ple of these are the Lorentz- 



Berthelot rules 






•A.B' 



^ ^ AA^ BB^ 



and 






AA 



BB^ 



(13) 



[8,9,1o], 



vjhich is based on London’s theory of dispersion energy 
•These rules give a better account of the second virial coeffi- 
cient of simple gases than the Lex'Jis-Rando.ll rule, yet there is 
real doubt ’whether any rule will yield accurate values of B^^ 

even under the best conditions of small, non-polar and spherically 
symmetric molecules. To make matters vrorse, very few experimen- 
tal measurements of second virial coefficients have been made on 
mixtures of simple molecules. Newer and more accurate functions 
such as the Kihara potential are presently coming under more and 
more study |^ 11 , 12 , 13 , l^j . 

V/hen consideration is on mixtures other than dilute gases, 
no exact relation is available to correlate the interaction ener- 
gies, Present practice |^6j is to extend the principle by 



considering an average interaction energy 



<UA 3 (r)> / C(x) = f[r/(7(x)] (l4) 

V/hile the virial equation provides the best theoretical basis for 
estimating the PVT properties of mixtures, it cannot deal satis- 
factorily V7ith any but low pressure gas mixtures. The inability 
to estimate third and higher virial coefficients generates an 
overall loss of accuracy. 

The thermodynamic properties of enthalpy and entropy are de- 
rivable from the virial equation |^4 , 15 f l6 , 17 J . From the First 
and Second Lavjs for a single-component, closed system dF=TdS - PdV 
and HmS+PV vrhence dH=;TdS+VdP, If this last is divided by dP and 
T is held constant, the relation ( ^ H/ 5 P)<j>=V-T( ^V/ ^T)p results. 
If Z=PV/RT, then ( ^ H/ 0 P)x= ( “^T^/p) ( 5 Z/5 T )p or after integration 

or 

(H 2 -Hp)^=-RT^^^( dZ/5 T)pd(ln P) (15) 



If P-,=0 and T =T/T^ , Pr=P/P. , then 

JL C C 

o rPr 

(H-H°)/RT^=-t2 J^{ ^Z/ ^T^)p^d(lnP^) 



(16) 



V7here (H-H*^) is the enthalpy change from P=0 to P=P at constant 
terraperature . Then for mixtures [ 3 j t Fqs. (2) & (15) give 

(Hj^-H°)/RT= 1/V [ Bm-TdB^dT] + l/V^ [ Cj,^-T/2(dCydT) ] + ••• (1?) 

From the Kaxwell relation (c)S/c)P)rp = (dV/5T)p and Z=PV/RT we 
have ( d S/5 P)qi=-R/P ^ Z+T ( 5 Z/ ^ T ) p j . Since the entropy blovrs 
up as pressure goes to zero, it is convenient to use the artifice 
of integrating this relation as a real gas to a pressure near 



zero where properties approximate those of a perfect gas. The 
equation is then integrated using perfect gas properties up to a 
convenient reference pressure, P , so that lov; pressure effects 

o 

on the entropy mutually cancel 1 4 . Then 



(S-S3 




(l-Z)d(lnP)+ln(P /P)- 

s 




d(lnP) 



(18) 



or in reduced terms 



(S~Sg)/R^^fl-Z) d(lnP^)+ln(P^^/P^) - T^^^( ^Z/5T^)p^<i(lnP^) 

(19) 

vjhere (S-Sg)/R Is the entropy change from 

T,Pg to T,P or to Then, Eq . (2) & Eq. (18) give 

/2V^+cy2V^ - T/2V^cydT] (20) 



Values can be determined for dB/dT by differentiating curves of 
B*(T) as a function of T" , using the relation dB /dT =(C/k/bQ) 



dB/dT . 

Finally to determine H and 

mm 

0 

and 8 ^ 3 = E - R 



the relations 



E 



•Inx^ + 



RlnPs 



( 21 ) 



can be used, vrhere is the enthalpy of component 1 at P=0 and 
S. is the entropy of component _i at reference pressure, P . 

The virlal equation is not the most convenient method of deter- 
mining the thermodynamic properties of 'mixtures unless computer 
facilities are available. If a virial equation has been formu- 
lated for PVT properties, however, it would not require much 
additional effort to program the derivatives of the equation. 



I 
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Pseudocritical Properties 

It vras stated earlier that the molecular energies of inter- 
action could be described by their dependencies on a universal 
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function as expressed by Eq , (l4). The molecular theories of 
Prigosine et al in vrhich explicit relations for the parame- 

ters £ (x) and (j (x) are developed in terras of their pure and 
mixture forns, ie, ^ aA’ ^ab’ ^AA’ ^AB’ 

Lennard-Jones or other ra-n potential forms, are formally equival- 



ent to the pseudocritical methods of Kay et al [sj. The relations 

kT/ £ {x) and V/(j'^{x) are equivalent to universal functions vrhich 

may be replaced by T/T (x) and V/V (x), vrhere T (x) and V (x) 

c c c c 

are pseudocritical properties of the mixture [ 3 ]. The pseudo- 
critical concept is based on the existence of a mythical pure 
fluid at the terraperature and pressure of the mixture vrhich exhi- 
bits the same properties as the mixture. The critical properties 
of this fluid are the pseudocritical properties of the mixtures, 

A sizable amount of literature has been vrritten on the deter- 
mination and use of the pseudocritical constants ^2,4,5,6,?j, 

Kay's simple mole fraction averaging rules |^2j are the most fre- 
quently used; in fact, the various rules for pseudocritical temp- 
erature all reduce under various simplifications to that of Kay, 

If T ET (x), Kay's rules are given by 
cm c 



•on.= S 



x.T . 

1 ci 
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and 



cm 






(23) 



1 Cl 



For T^j^, Kay’s rule produces values vrithin 2% of other rules for 
0. T^^/T^gC 2.0 and 0 • 5 2 [sj. For values of P^ dif- 
fering by 20,^ or more, the modified Prausnitz and Gunn rule |^9j 
should be used. 



P =Z RT /V =R 
cm cm cm cm 
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Kay*s rule has been verified to predict the true critical temp- 
erature of N^-Og mixtures to within 10°K, 

Using the pseudocri tical rules for mixtures of nitrogen and 
oxygen, the author computed enthalpy deviations by the method sug- 
gested by Reid & Shervrood. A summary of this method f ollov:s : 

1. Pseudocritical constants for the mixture were calculated from 
the follovjing pure component values; 



T 

c 

P 

c 

V 

c 

Z 

c 



Nitrogen 




Oxygen 


126.26°k 


[ii.i'f] 


15^.77°K [12,13] 


33,5^ atm 


[ll.l^] 


50.14 atm [12,13] 


90.07 cm^, 


/mole [15] 


.075 cm^/mole [ 12] 


.291 [s] 




.292 [s] 



2 . 



3. 



Ba.se enthalpies, ^'or nitrogen and oxygen were selected 

based on the enthalpies of the perfect gas at 1 atm at the 
normal boiling point and triple point respectively |^ll,12j , 
The ideal state enthalpies at the various isotherms of inter- 
est, H°, are then calculated from 




(25) 



where . 121 j/mole°K and 29,359 j/mole°K [l6]. 
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4 . 



At a oarticular temperature the perfect gas mixture enthal- 

pies, H°, may be calculated from 
m 
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5. 



The mixture entVialpy at the system pressure is then 



( 27 ) 






m 



m 



vrhere the isothermal enthalpy deviation can be deter- 

mined from suitable tables ^lyj , 

Reid and Sherv:ood predict errors of lOt in the vapor phase and 
20% in the liquid for this method [sj. This author carried out 
calculations at 84°K and 1 atm and at 112°K and 10 atm. In the. 
vapor phase for the case of 84°K and 1 atm, although the calcula- 
ted values differed from those obtained using a mole fraction 
average of the pure component enthalpies by as much as 35 j/mole, 
the results vjere linear and had the same slope (Pig, 4a). In the 
liquid region at this temperature the results shovred a small neg- 
ative deviation from linearity (experimental evidence indicates a 
positive excess heat) vrith the absolute difference increasing 
v;ith increasing ^mole O 2 (Fig. 4b). At the higher pressure the 
calculated vapor isotherm was no longer linear, but exhibited a 
sizeable negative excess heat. Likei-rise the liquid Isotherm at 
this pressure displays a somevrhat larger contrary behavior than 
before. All isotherms are extrapolated into the wet vapor region. 
These calculations seemed to be undesirable for their deviations 
from known behavior, ie, negative excess heat. The tables offer, 
however, a relatively simple computational technique and are ap- 
plicable to the liquid phase as well as the vapor. Perhaps better 
values of the pseudocritical properties would produce more logi- 
cal results, but for the present, simple mole fraction averaging 
seems to provide more satisfactory answers. 

While tables such as those of Lyderson et al do not. exist for 
entropy deviations from the perfect gas base, a relationship can 
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Figure 4a. 
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Figure 4a & 4b. Liquid and vapor lines determined by variou 
methods. Lashed lines join pure component values’. 



be viritten betvreen properties of other values which are tabulated. 
Eq. (18) can be written 

(S - S 3 )/R=(H - H°)/RT - ln(f/Pg) (28) 

when P=Pg values of both (H - H°) and (f/Pg) have been tabu- 

lated. 
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Empirical Squatlons of State 

Perhaps the most underdeveloped of the techniques for the 
determination of the properties of mixtures is that of a workable 
equation of state, particularly in the tv;o-phase region. 



The Benedict-'.'/ebb-Rubin equation v;ith eight constants 

has been fruitfully tested, but unfortunately is applicable only 
to light hydrocarbon mixtures. 

Some effort has been made to determine coefficients for the 
two constant Redlich-Kx-Jong equation as follows: 



It is reported that some success has been obtained in using 
reduced temperatures and pressures derived from Kay's rules in 
the generalized Beattie-Br idgman equation 



Recent ef essful in treating 






(29) 



X'There for a mixture 




and 







(30) 




gas mixtures of CO^ and M 2 with mixture rules for the Martin-Hou 



equation of state. 

5 




= A. + B.T + C,exp(k'T/T^) (31) 

J. X J. c 
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None of the equations of state have been extensively tested with 
mixtures, v;ith the possible exception of the BV/R equation for 
light hydrocarbons, and all require machine handling. In the 
derivation of thermodynamic properties using Sqs, (15) & (18) 
the results will be only as good as the equation used, and some- 
times v:orse, because of the errors involved in the determination 
of the derivatives of the equation. 
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Partial Properties 

It vras mentioned earlier that use might be made in some cases 
of partial properties in determining mixture properties. Partial 
molar properties define the properties of components in a mixture 
as opposed to those of the pure components and represent a part 
of the total property of the particular mixture. In particular, 
a property is defined as 








(32) 



v;here R is any extensive property and n^^ is the number of moles 
of component _i in the mixture. These relations are governed by 



I 




the Gibbs-Duhem equation 



y~. dH^ - 



( 33 ) 



If R = R/n = x.R. , 
i 



then it can be sho 



[^] 



that 



R^ = R - Xg(dR/(^Xg) and Rg = R - x^{ dR/ (3^) 

The Lewis relation for fugacity as it is related to the 

chemical potential can be expressed for a inixture as 



where is some function of temperature only. As the pres- 

sure approaches zero, f approaches as opposed to the chemical 
potential which blows up as can be seen from ( d jti^y ^R)rj. ~ . 

At constant temperature and composition Sq, (35) can be written 




RT ln(f^/f^x^) 




Vj^)dP = RT In 



(36) 



vrhere Y i defined as the activity coefficient. For ideal sol- 
utions this reduces to' the Levris & Randall rule [2 j 

f^ = x^f^ such that = 1 . (37) 

While in an ideal solution it is only necessary that the inter- 
m.olecular forces be the same for all components, in a perfect gas 
the molecular forces must be zero. In this case the pure compon- 
ent fugacities are equal to the pressure and Raoult's law results 



I 
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f = P, = xYp 



i i i 



(38) 



It should he mentioned that the Levris -Randall rule generally 
provides a good approximation for most real non-polar gas mix- 
tures at any pressure as long as the component in question is sim- 
ilar to the other component or present in excess of the other com- 
ponent , 



If the fugacity of the pure components in the mixture can be 
determined, the enthalpy deviation representing the non-ideality 
of the mixture can be determined from the relation 



The mixture enthalpy then, follows from the partial property and 
perfect gas base state relation 



vrhere the pure component perfect gas values are determined as be- 
fore . 
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Excess Functions 

Before discussing deviations from ideality, it would do well 
to define somev;hat more fully the ideal mixture. One definition 
j^lj states that it is a mixture whose components have chemical 
potentials given by 



M^(P.T.x) = Mi(P.T) + RT Inx^ 



(^ 1 ) 



where ^(P,T) is the pure component potential at the pressure 
and temperature of the mixture. It follovrs that for an ideal mix 
ture 



G = 
H = 
S = 
V = 





Inx^ 




XiH, 



X S. R 

i i 



x.V. • 

1 1 



E 



x^lnx 



i 



(^^2) 



The log sums are referred to as the ideal free energy of mixing, 
G^, and the ideal entropy of mixing (Gibbs paradox), S^. The 
ideal free energy of mixing is alvrays negative and that of entro~ 
py always positive. The ide'al enthalpy and molar volume of mix- 
ing are equal to zero. 

Excess properties may be defined in these terms as the dif- 
ferences betvreen the values of the properties in the real mixture 
and the values of the properties given by Bqs . (^1-2) for an ideal 
mixture of the same composition. 



2 ? 



It has been seen that the chemical potential may be related 
to the activity coefficient, so that -for a real mixture at low 
pressures 

/^.(P,T,x) -//.^(P.T) = RT (43) 



The partial molar properties may then be given by 

- RT^(dln/^?5T)p 

- R - RT<d lnXj^/^T)p 

+ RT(^ln7./^P),p 



(44) 



The excess thermodynamic properties may then be obtained by sub- 
tracting* the ideal quantities from the partial molar quantities. 




V® = RT Zx.(e)ln/./c)P)T 
The excess functions are related by 
- TS® = + T(dG^/5T)p 



(45) 






(46) 



Considerable effort has been m.ade in the past decade in the mea- 
surement of excess functions. In particular, regarding ^£"^2 
mixtures, Knobler et al ^2,3»4j measured H"" and V® at 7?°K; Wil- 
helm & Schneider [s] » at 77.5°K; Pool et al [^j, g"" at 83.8°K; 
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and Duncan & Staveley r at 63 . 1 ^+°K. In addition, values 
E 

of G have been calculated by some of these authors fx'om the 



[a]. 



Armstronp: et 



vapor pressure measurements of Dodge & Dunbar 
al ^9] » find Din [l^j , These results are plotted in Fig, 3 ?- for 
equimolar liquid mixtures. 

It can be shown [^ij that the variation of vapor pressure at 
constant temperature can be related to the excess Gibbs free en- 
ergy, Vapor pressure and composition can be measured in an 
equilibrium still, a large number of which have been devised 
^ll- 24 j , The principal problem in all stills is the difficulty 
in eliminating vapor condensation before it reaches the vapor 
sampling point and of preventing liquid droplets from being car- 
ried over to the same point. These problems can be eliminated by 
methods vrhich utilize only the total static vapor pressure as a 
function of liquid composition. Then the vapor composition may 
be found analytically using the Gibbs-Duhem equation. Several 
methods have been used to measure the total vapor pressure [25- 

29]. 



If isothermal vapor pressure measurements are knovm, values 

E 1 V 

of G may be derived from knowledge of P, x , x and the second 

virial coefficient. If the easier experimental procedure has been 

followed, vrhich measures x^ and not x^, there are several graphic, 

algebraic, and machine methods vzhlch may be used to solve for G 

1^27,30-35]. These solutions are more complicated than those used 

when both compositions have been measured. Other methods have 

been devised for the calculation of the liquid composition, if 

only the vapor composition has been measured ^ 36 - 3 b] . Finally, if 

all four properties are knovm, then the thermodynamic consistency 
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of the results may be tested v.'ith the Duhem-Kargules equation j^39 
The most vridely used test for consistency is that developed 
by Coulson & Herington [4o] and by Redlich & Kister [^l] , If the 
volatility ratio is defined by 0( = Xg/xJ / x^/x^ where B is the 
more volatile component, then the difference in excess chemical 
potentials of tvro components can be given as 






B 




HT ln(7g/7;,^)l 



= RT ln(CX P^/Pg)+(P-P^) 



<''1 - V> - - PbHvJ - Bgj,) + 2P(6B)^g(xf- xf) (47) 

By intGgrs-ting this BxprGssion ovGr X-nt where is zero 8.t x=^:0 
and x=l, V7e obtain 



flniCK PyPgjdxl 



= 1/Rt| [(P - PgXvl - B^g) 

-2P(6B)^3(x.f-xf)]dxl 






(48) 



If the results are consistent, the volatility ratio will equal 
the pure component pressure ratio, ie , 0( =Pg/PA* the composi- 
tlon V7here G is a maximum or a minimum value in a non-ideal mix- 
ture in equilibrium with a perfect gas. Therefore, by plotting 
^ function of Xg the areas above and below the axis 
v/111 cancel if the results are consistent, 

« 

Unf ortunately , it has been demonstrated [ 28 ] that this test 
fails for liquid mixtures in the vicinity of their normal boiling 
point Sj because it assumes tnat the vapor behaves as a 'ooi*rect 
gas, Rowlinson |^lj modified the method by noting that ((3 b) 

A.B 

made no significant contribution and. that P could be replaced by 
Raoult's law pressure, x^P^ + XgPg, so that 
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^ ln(0( P,Y/Pg)dX3 - [(P^ - P 3 )(v| - vj - - Bgg)J /2RT 

an explicit relation irhich is convenient to use. 



(49) 



The most widely noted method of determining G® from measured 
isothermal values of total pressure and liquid composition is that 



of Earke 



:r [41]. 



This author has made calculations of for 
Ng-Og mixtures at intervals of 10°K from 80°K to 120°K using 
this method. The procedure is outlined for interest in Appendix B, 

Rowlinson [l] has indicated the difficulty in defining the 
excess free energy at vapor pressures greater than 2 to 3 atm, 
because of compressibility effects. This can be seen in Fig, 32 
where the values of G® calculated by the author have been plotted 
with those calculated by other investigators at lower tempera- 
tures. It should be noted, however, that the slope of the values 
betv;een 100 and 120°K (dashed line) is identical to that at lo’wer 
temperatures (solid line), and one wonders if perhaps a method 
for calculating G^ including the third virial coefficient might 
not displace these values into line. 

At any rate, studies in the vicinity of the two-phase region 
at pressures higher than atmospheric are virtually non-existant 
[ 1 , 43 , 44 ]. Even the molecular theorists have contented themselves 
vrith predictions of excess functions in the liquid phase at neg- 
ligible pressures or when venturing into the field have confined 
themselves to predictions of relations for the critical constants. 
In the gas phase preoccupation has been with determining an equa- 
tion of state [45] and almost no experimental vrork has been done 
to determine excess properties of gases. Only very recently have 
badly needed measurements on dense gases been attempted [46,4?]. 
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It had been generally felt that the molecular theory 

vjould stand up vjell in the estimation of excess functions of mix- 
tures of condensed gases, that is, liquid mixtures of simple mol- 
ecules at temperatures near the boiling points of the mixtures. 

In the past decade, however, experimental vjorh by Beenakker et al 
[ 2 , 3 .^] and Staveley et al [ 6 , 7 ] have shown serious inconsisten- 
cies and plaguing weaknesses in the molecular methods. They have 
•found that values of excess energy calculated by the molecular 
theory of Prigogine et al [^^ 3 ] are strongly dependent on both 
the set of potential param.eters chosen and on the pixre component 
selected as a "reference" substance. 

The most probable^ cause of failure of the theory is consi- 
dered to be the Lennard- Jones pote^;itial [ 6,7i^S,49,5oj particular- 
ly at the close molecular distances involved in liquids. Still 
another cause may lie in the innaccuracy of the interaction par- 
ameter combining rules. At any rate, the spectacular differences 
in the value of excess functions calculated from molecular theory 
from those measured experimentally is enough to preclude at this time 
consideration of values formulated by -this means. 
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Calorimetry 

The measui'ement of the excess .heat is of such importance as 
to deserve a separate section revie’.:ing devices and methods pres- 
ently available for its determination. The two principal instru- 
ments are the adiabatic and Isothermal calorimeters. The adia- 
batic calorimeter reflects pos itive. excess heat, if the tempera-- 
ture drops on mixing two pure components in an adiabatic contain- 
er. If on mixing or afterwards heat is added to cancel out the 
termperature drop, the am.ount of excess heat can be measured. 

Heat added during mixing reduces corrections necessary for heat 
transfer to the environment, A similar determination may be made 
for negative excess heat with its indicative temperature rise. 

The adiabatic calorimeter has been bothered by a variety of 
problems including imperfect heat isolation, incomplete mixing, 
inaccurate measurement devices, and spaces in which vapor can 
form, The last problem is complicated by the need to allow for 
volume changes on mixing. 

Isothermal calorimeters have the advantage of not requiring 
corrections for heat losses or gains and are able to measure a 
greater temperature range accurately. 

Many such devices are mentioned in the literature and some 

* 

of the more important are listed here. Some early calorimetry 

*o * 

at temperatures of about 100 K was done by Howlinson 5: Malcolm 
^ij, A versatile instrument vihich can be used either isothermally 
or adiabatically has been developed by Nicholson . A glass, 
highly accurate, apparatus is described by Larkin 5: McGlashan 
vrhich has eliminated vapor spaces, provided complete mixing and 
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allov:ed for expansion and contraction through a compensating 
piston, Lilcevrise, Mrazek & Van ITess in their equipment con- 

trol the volume with close clearance pistons. Another novel con- 
trol on volume change is the metal aneroid diaphragms of the de- 
vice of Brovm & Pock • Work has been done vrith Bunsen iso- 
thermal calorimeters by some investigators [6,?], Two recent 
calorimeters were specifically designed for measurements of con- 
densed gas mixtures 
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Latent Heat 

The latent heats of vaporization for mixtures are covered by a 
multitude of definitions The most obvious are the heats 

required for vaporization or condensation at constant pressure, 
constant temperature or constant composition. 

One method of calculating the latent heat is by employing 
variations of the Clausius-Clapeyron relation 1,3-8 . One of 
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these variations expresses the differential Vieat of vaporization 



/\v^^dP - (/\K^^/T)dT - (x^ - x^)( x^^)p^^dx^ = 0 (50) 

for a small amount of vapor phase out of or into a large amount 
of liquid. Similarly, the differential heat of condensation is 

A^ivdP - (AH3^^/T)dT - (x^ - x^)( d2o^/6x^2)^^_^dx^ = 0 (51) 

for a small amount of liquid phase out of or into a large amount 
of vapor. The volume change is given by 

Av^P = - (x^ - x^)( ^)V^/5x^)p^^ (52) 

and the heats are 

Ah^^ = - (x^ - x^)(dnV^x^)p ^ (53) 

and represent the heat required at constant temperature and pres- 
sure to obtain one mole of vapor phase from a large quantity of 
liquid or to condense one mole of vapor into a large amount of 
liquid. Similar expressions apply for A A^P^* -ven tfith 

this seemingly straightforv:ard approach there is no simple method 
of solution except with assumptions vrhich ignore the effects of 
non ideality (see Graphical Kethods). 

Corresponding states theory yields an imposing array of esti- 
mation techniques for pure component latent heats |^7j . Calcula- 
tions were made by one of these methods |^9j, alleged to be among 
the most accurate and simple to apply, by adapting it to mixtures 
by use of the pseudocrl tlcal constants. This em.pirically derived 
relation can be expressed 
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<A«vl>b ' [Tb<5 - 2.17)] /(0.93 - (5'0 

vfhere is the nornal hoilin.<^ point and is the isothertral 

differential heat of vaporization. The results, calculated at 
one atm, were less than satisfactory as can be seen in Fig. 31 • 

In fact, because of the large latent heats as compared to the ex- 
cess heats Involved, it is a more accurate estimate simply to use 
the mole fraction average of the pure component latent heats, 
particularly at low pressures. 

One experimental determination of the latent heats of N 2 -O 2 
mixtures exist j^ioj » ^ notable number vihen so fevr for any binary 
mixture can be found, Dana's measurements at atmospheric pressure 
are plotted in Fig. 31. They represent the integral heat of vap- 
orization at constant pressure, that is, the heat required to 
change one mole of liquid to the vapor state under conditions of 
constant pressure and composition. 

No matter hov: the heat of vaporization is determined, hov7- 
ever, it is only as useful as the isotherm or phase boundary from 
which it will be based. 

References for Latent Heat 

1. Gamblll, Chem. Sng. , 65 , 13?, 1958. 

2. Fischer, Ann, der physlk. , 1?, 209, 1933. 

3 . Dodge, Chem, Eng. Thermo*., McGraw h!lll, NY, 19^^. 

4. Edraister, AIChB J. , 1, 38 , 1955. 

5 . Hobson & Weber, Petrol. Process,, 12, 15, 1957. 

6. Bosnjakovic, Blackshear (trans). Tech. Thermo. , Holt, Rine- 
hart 8c Winston, 19^5 . 
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7. Reid <Sr Sherwood, The Props, of Gases & Liquids , McGrav; Hill, 
NV, 1966 . 

8. Ruheirann, The Separation of Cares , Oxford, 19^9* 

9 . Reidel, Chem. Ing. Tech., 26, 6?9» 195^. 

10. Dana, Proc. Am. Acad, Arts & Scl, 60 , 24l, 1925. 

Graphical Methods 

« 

The earliest graphical techniques for determining the thermo- 
dynamic properties of Ng-O^ mixtures v;ere devised at the Kamer- 
lingh Onnes Laboratory at the University of Leiden Keesom 

& Tuyn developed a general construction for H,S,x surfaces of the 
mixture in the tv;o-phase region at atmospheric pressure. The 
property of the mixture that is the key to the geometric method 
is the transverse fold in the H ,S , x-surface representing the vret 
vapor region. A tangent plane across the fold contacts the points 
of equilibrium or coexistant phases. These points were called 

connodes ; the isotherm joining them, the nodal line; and the 

# 

locus of these points, the connodal curves (vapor and liquid 
curves ) . 

Basically, their procedure was to locate known points, such 
as the boiling point of oxygen on an H-S diagram, and with the 
aid of pure component heats of vaporization at I'atm, the defi- 
nition of enthalpy and entropy in terms* of specific heat, ideal 
mixture rules, and Dodge and Dunbar equilibrium data [^3j. plot 
the vapor line which is then easily projected into the S-x plane. 

Using the previously mentioned principle of tangent planes, 
the liquid line was determined in the H-S and S-x planes by 
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graphically constructing the line representing the plane tangent 
to the connodal points, assisted by Dana's heats of vaporization 
at constant composition [^] • The projection in the H-x plane v;as 
then easily established. 

The procedure then recommends extrapolation in the H-S plane 
of the lines of constant composition to higher pressure using 
Second Law relationships and the virial equation of state, ie. 



dH = V - T( d V/5T)p dP and 

r f X 

dS = V/^T)p^^dP 



M55) 



which become v:ith the application of the equation of state 



H = (P 2 - ^2^^® “ 

S = R In P]^/^2 *■ ^^2 " ^1^ dB/dT 



(56) 



The ideal mixing rules locate the isocomposition lines between 
those of the pure components, and the equilibrium data, fixes the 
end points at the vapor line. The form of the liquid line is 
then taken to be the same as that for .one atm. 

Lerberghe proposed a variation on this method of his 

colleagues which requires only the use of the T-x equilibrium 
diagrams. From thermodynamics and analystical geometry, expres- 
sions were developed for the enthalpy and entropy difference be- 
tween the connodal points. 



= (N + KT)(x"^ - x^) and 
= (L + M) (x^ - x^) 



M57) 
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vrhere 



V V fy V . . V -| 

vrhere L = - S 3 + R In I ( 1 - x )/x 

M = -Rt[i/(1-x^) + l/(x^)] dx^/dT 



(58) 




where the A and B are the pure component properties at the pres- 
sure and temperature of the mixture, extrapolated as necessary 
into the v;et vapor region. 

Thus by fixing the vapor line as was done by Keesom and Tuyn, 
Lerberghe’s method locates the liquid line with the help of the 
equilibrium data. Note that the requirement exists to perform 

V 

the difficult graphical- tangent determination for dx /dT or to 
derive an equation for the vapor line on the T-x diagram by curve 
fitting and then taking its derivative analytically, Lerberghe's 
graphical measurements are reflected in Fig. 31 For the liquid 
line at one atm. 

Several years later Bosnjakovic ^ 6 j nicely summarized vari- 
ous graphical proce.diires . His method suggests that in most cases 
the excess functions in the vapor phase can be neglected, at 
least in moderate pressures, and the isotherms in the H-x and S-x 
diagrams may be drawn according to the ideal mixture rules. Then, 
if the temperature-vapor composition relationship is known from 
equilibrium data, the vapor line is also determined in both 
planes , 

In the liquid region, if the excess properties are knovm, 
then the isotherms m*ay be drax-rn and the liquid line again deter- 
mined as above for the vapor line. If the excess heat is not 
knox^n, some estimation of the latent heat of vaporization must 



be made. The Clausius-Clapeyron relation (Eq.-51) makes a good 
starting point. From the definition of Gibbs free energy the 
second derivative is given by 

T( (59) 

Then for an ideal vapor mixture 

(6^G^/^x^2) RT/x^(1-x^) (60) 

r , i 

Measurement of the slope of the T-x diagram graphically or analy- 
tically yields (^x^/^T)p, Enough information is then' available 
to solve Eq. (51) for A constant pressure and temperature. 
The vapor isotherm is then extended to intersect x^ and A ^ 

Iv 

is plotted from that point (Fig. 5). Repeating this process 
V7ith isotherms of the vret vapor region yields the liquid line 
(Fig. 31) • Obviously the method is sensitive to the behavior of 
the two differential terms of Sq. (5D* 




In the case of the S-x diagram, .after the Isotherms have been 
drawn according to the ideal mixing rules and the vapor line de- 
termined from equilibrium data, the liquid line may be determined 
from the relation 



Ah^ /T = 
iv 







P,T 



(61) 



The value of /\h^.^/T can be obtained from the H-x diagram and 
plotted from the intersection of the tangent to the vapor line at 
x^ with the equilibrium liquid composition, x^ (Fig. 6). 

A secondary method of constructing the S-x diagram when the 
H-x diagram is available utilizes an integration procedure, which, 
needless to say, is less prone to graphical error than the dif- 
ferentiation technique. This method uses the relation that, for 
a reversible vaporization process at constant pressure and con- 
stant liquid composition. 



(s'" - S^)^l 






= J 1/T (dH)^l 



( 62 ) 



The right hand side of the equation may be easily integrated 
graphically using any of the several numerical integration rules 
and the difference measured from the known value of S^. 



References for the Graphical Methods 
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Rinehax't & V/lnston, 1965. 

7. Ruhemann, The Semratlon of Gases , Oxford, 19^9. 

The Thermodynamic Properties 

The thermodynamic properties of binary mixtures of nitrogen 
and oxygen considered to be of the most value for use in engineer- 
ing calculations are those of enthalpy and entropy, A useful 
format for the presentation of these properties is the graphical 
display, in particular the enthalpy-composition and entropy- 
composition diagrams at constant pressure. Accordingly Figs, 

7-30 have been prepared for mixtures at pressures of 1, 5. 10, 
and 20 atm. Enthalpy and entropy have been represented as func- 
tions of composition in the liquid line region, the vapor line 
region, and the vret vapor region. 

Values for the pure component oxygen v;ere taken from Stev:- 
art's machine tabulation ^ij based on a complex empirical expan- 
sion, Those for nitrogen came from the work of Strobridge ^2j 
whose machine tabulation is likewise based on an empirical equa- 
tion of state. 

The molecular v;eight of nitrogen v;as taken as 28,016 
grams/mole ^2j and one calorie as 4.184 abs, joules, by which 
conversion factors the tabulated values of Strobridge v:ere con- 
verted to abs, joules/nole. 

For the sake of graphical convenience the properties of 

o 

nitrogen and oxygen were equalized at 91 K at 1 atm, creating an 
oxygen enthalpy scale factor of 4245.3 j/mole and a nitrogen en- 
tropy scale factor of 81,12 j/mole. The values of entropy are 



relative only, since the nitrogen scale is not. based on zero en- 
tropy at absolute . zero as is that of oxygen. Property values of 
the pure components within the tv;o phase region have been linear- 
ly extrapolated. 

Values of the entropy of mixing (Gibbs paradox) are derived 
from the relation 



,m 



R 




Inx^ 



( 63 ) 



—2 

where R is the universal gas constant equal to 8,2057^ x 10 
L-atm/mole-°K j^lj , 

Isotherms in the superheated vapor region vxere derived as 
mole fraction averages of the values for the pure components, ie, 
from the ideal mixture rules. The excess heat and entropy func- 
tions for the liquid isotherms vrere derived by the author from 
the data of several investigators (see Discussion of Results). 

The vapor and liquid lines vrere plotted using the isotherms and 
the temperature-composition data of Dodge and Dunbar 3 • 



References for Thermodynamic Properties 

1. Stevrart, Ph.D, Thesis, Univ. of Iowa, 1966. 

2. Strobridge, NBS Tech Note 129, 1962. 

3. Dodge & Dunbar, J, Am, Chem. Soc., 49, 591, 192?. 
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DISCUSSION OP RESULTS 



The choice of the method of determination and the form of 
presentation of the results v^^as arrived at after consideration of 
all the previously discussed techniques. No attempt was m.ade to 
formulate a virial expansion or other empirical equation of 
state, because of a general lac'k of the PVT data required for 
such a formulation, A corresponding states theory estimation 
technique, using pseudocritical constants and the enthalpy devi- 
ation tables of Lydersen., was conducted, with the results already 
mentioned. The negative excess heats calculated by this method 
are thermodynamically inconsistent, and this procedure was re- 
jected. Using pseudocritical constants, a modified relation for 
the latent heat was calculated, which when plotted produced the 
least agreeable of all liquid lines determined (Fig, 3D* The 
purely graphical approach of Lerberghe produced better results, 
except at high molar concentrations, v^hlch plotted extremely 
erratically,. Both a mole fraction average of the heats of vap- 
orization of the pure components and a mole fract.ion average of 
the pure component boiling points produced better results in 
that order, than did the complex estimation techniques mentioned 
above, A plot of the results obtained by these methods at 1 atm 

o 

'together v;ith the method finally chosen is given in Fig. 31 • 

I 

In the end it was decided that the calculation of excess 

^ . • 

functions provided the best determination of the liquid line, and 
therefore of the isotherms in the liquid region. The liquid 
line determined by this method shows very good correlation with 
the experimental evidence of Dana, as can be seen in Fig. 31. 
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0 0 — Rei^el pseudocrit. const, determination of(AHvl)b. 

-A A — A-Lebirghe graphical determination of (AHvI). 

Mole fraction avg. of pure component liquid points. 

Mole fraction av^. of pure component (A^x). 

■ Excess heat determination from this work. 

0000 Experimental data of Dana. 

Figure 31* Liquid line for N2“02 mixtures at atmospheric pressure. 



A presentation of the excess Gibbs free energy as calcula- 
ted by previous investigators is shown in Fig, 32. The excess 
Gibbs free energy is related to the other excess functions by 

= G^ + TS^ = G^ - T( ^>G"/^T)p (46) 

p 

It is clear from Fig. 32 that the slope of G with respect to 
termperature is very linear in the lovr pressure region and has 
been measured as 0.08 cal/raole °K. With this slope the excess 
heat of the liquid v;as calculated to be l6.0 cal/mole at equimolar 
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composition and lovr pressure. Only one experimental calculation 
of excess heat for ^2**^2 fixtures has ever been made [ij. 

These measurements by Knobler et al of the excess heat at 77°K 
are plotted as Fig, 33. These results can be expressed by fhe 
fitted polynomial approximation |^2j . 



= %2 ""02 



[42 + 10(2xq^ - 1) + Z0(2Xq^ - l)^]cal/mol 



O2 



(64) 



which yields an excess heat of 10,5 cal/mole at x=0,5. A concur- 

p 

rent calculation of O'" was made and can be given by 

= Xj ^2 V 2 + ^(2xq2 - 1) + ^(2Xq^ - D^jcal/mole. ( 65 ) 

This gives an excess Gibbs free energy of 10,0 cal/mole at x=0,5 

or a value nearly the same as that of the excess enthalpy. Since 

. E EE 

the temperature derivative of G is negative, then H :=> G , and 

from the vxeight of evidence by considerably more than 0,5 cal/mole, 

■It v:as decided that l6,0 cal/mole was a reasonable value for the 

excess heat at equimolar composition In keeping with Fig. 32, 

t 

and the good fit obtained in plotting the liquid line; Knobler *s 
data was then scaled up to this result. It should be noted that 
in agreement vjlth prediction [^ 3 ] the excess heat is asymmetric 
'with respect to the comiposit ion , while the curve of excess free 
energy is symmetrical. Values of the excess entropy were then 
calculated from Eq. (46), 

When the calculated values of excess heat were applied to 
isotherms in the liquid region at one atmosphere, the resulting 
liquid line showed excellent agreement with the experimental 
data of Dana’, It was decided, therefore, that this technique 
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Figure ^2, Excess Gibbs free energy as a function of temperature 
for equimolar liquid mixtures of N 2 -O 2 . See Appendix C for data' 
sources. 
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Figure 33. Excess heat and excess Gibbs free energy as a function 

•p 

of composition from the data of Knobler et al. H fitted to re- 
vised equimolar avlue of this thesis. 
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at atmospheric pressure was the best available. 

The Integration method of Eosnjakovic was used to calculate 
(S - S )^ 1 . The liquid line thus determined for the S-x dia- 
gram was in excellent agreement with that calculated from values 
of the excess entropy. 

In the vapor phase it had been decided that most authors 
generally agreed that there was negligible heat involved with ' 
mixing at low pressures including atmospheric, so that the best 
approximation of isotherms was a simple mole fraction average 
corresponding to the mixing rules for an ideal mixture; this in- 
cludes the entropy of mixing for the S-x diagram vapor line. 

The major problem then, concerned pressures greater than 
atmospheric, that strangely Ignored and abandoned region vrhere 
investigators have feared to tread. It would be foolhardy and 
certainly not justifiable to make unwarranted assumptions about 
conditions at higher temperatures v;here pressures exceed- several 
atmospheres, where the third virial coefficient, the compressi- 
bility of the liquid, and the volume change on mixing are all 
important factors. 

For the best available presentation it was decided that the 
excess functions of heat and entropy already calculated would be 
used to determine the isotherms of the liquid and mole fraction 
averaging the vapor isotherms at pressures and temperatures be- 
low the critical. This decision is not wholly v;ithout justifica- 
tion. It is not unreasonable to assume that the excess thermo- 
dynamic properties of the liquid are relatively independent of 
pressure [4], principally because of the relatively small volume 
changes in the liquid state. Additionally, Fig. 32 gives some 
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small hint that deviations at the highej? temperatures associated 
with the higher vapor pressures are not so great from the more 
accurately calculated linear results at lovxer temperatures. 

The incompressibility reasoning does not apply so well in 
the dense gas region, but here too it is felt that the ideality 
of mixture assumption is valid enough vxhen compared to other 
available data. The values of several investigators for air 

were compared at the pressures in question at several 
isotherms near the dew point and compared to the 79.^wole N 2 mix- 
ture of nitrogen and oxygen isotherms calculated by mole fraction 
averaging. Although the deviations varied widely, a positive 
excess heat with increasing pressure might seem to be indicated, 
the effects of argon in the mixture and the particular investi- 
gations notwithstanding. The deviations are not large enough nor 
the.ir pattern clear enough to justify guesswork at possible val- 
ues of excess heat. 

References for the Discussion of Results 
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1959. 

4, Bosnjakovlc, Blackshear (trans.). Holt, Rinehart & Winston, 
NY, 1965 . 

5 , Michels, Vfassenaar & Wolkers, Appl. Sd. Res., A5, 121, 1954. 

6, Lansbaum et al, AIChB J., 1, 302, 1955. 

7, Hilsenrath et al , NBS Circ. 564, 1955. 
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CONCLUSIONS 

The intention of this thesis v;as to survey the literature 
with regard to the thermodynamic properties of binary mixtures of 
nitrogen and oxygen in the two-phase region and to analyze the 
resulting data. From^the sizeable collection of available print- 
ed matter which concerns itself directly and indirectly with the 
subject, it would seem an easy matter to make an analytical de- 
termination of the common thermodynamic properties of such rela- 
tively simple mixtures, A wide variety of estimation techniques 
do, in fact, exist, and while those have been presented here, 
most, if ■ not all, of them are either unsatisfactory or at least 
only partially useful. 

Specifically, the virlal expansion and other equations of 
state with their associated derivatives would be most satisfac- 
tory if experimental PVT data existed with which to determine 
their coef f j.cients , Nolecular theory, similarly! has great po- 
tential, particularly for liquids and -dense gases, but suffers 
from the uncertainty involved in its. basic parameters. The tech- 
niques involving pseudocritlcal constants give good rough esti- 
s 

I mates, but do not provide the accuracy required for even engi- 

t 

* neering calculations. The excess func,tlons provide a rapid and 
relatively easy method for determining the thermodynamic proper- 
ties, Presently available data on these give good results at 
pressures near atmospheric. 

All of the aforementioned methods have one thing in common. 
None of them have been developed for use in the high pressure 



area of the tvxo-phase region, VIhere no one has gone before, hoX'X- 
ever, any determination, no matter how simplified, idealized, or 
approximated it may be, is better than none at all. Therefore it 
was decided that by taking the excess functions independent of 
pressure in the liquid region and using ideal mixture rules in the 
vapor phase, satisfactory thermodynamic diagrams could be con- 
structed , 



RECOMMENDATIONS 

Under this heading of Recommendations should be listed the 
solutions to the problems posed by the Results, At the risk of 
redundancy and speaking only for the mixtiires in question these 
inadequacies are repeated here, 

A single equation of state for all mixtures vrould be the 
most satisfactory solution as it vxould provide in one analytical 
expression and its derivatives the measureable and derived prop- 
erties for the infinite number of mixtures, A .solution of this 
type would require an extensive experimental program to deter- 
mine a sufficient number of PVT data to permit the evaluation 
of appropriate coefficients. 

Another satisfactory program would be to make direct measure- 
ments of the excess heat function using one of the calorimeters 
mentioned previously. This is the only satisfactory means of 
determining this property and at high pressure the only means. 

Then since the vapor pressure method of determining G is only 
satisfactory at essentially atm.ospherlc pressure, the excess en- 
tropy function could be determined by the integration procedure 
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discussed earlier under Graphical Methods, A program of this type 
V7ould not be all inclusive in its results and would necessitate 
taking measurements at sufficiently close intervals of temperature 
to permit reasonable interpolation of tabulated data or accurate 
drafting of thermodynamic diagrams. 
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Appendix A 

1, Temperature-composition diagram, 1 atm, 

2, Temperature-composition diagram, 5 atm, 

3, Temperature-composition diagram, 10 atm, 

4, Temperature-composition diagram, 20 atm. 
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Figure 54 




Figure 35 
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Figure 36 
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Figure 37 
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APPENDIX B 

Barker’s method: 

The most widely noted method of determining G from measured iso- 
thermal values of total vapor pressure and liquid composition is 
that set out by Barker, 

The excess chemical potential is given by Eq, (47) which can 
be written 

= RT In ^,2 = RT ln(P)c;;2/’'N2='J2’ ‘''n2 ' ‘ 

*■ f’6N202<’'02>^ 

II E 

with a similar equation for Mo2* 

It follovrs that 



^ - If K2%2 ’’’ f 02^02 



(67) 



where Pj^^ ^03 ^®P^®^ent partial pressures if the liquid is an 
•ideal solution but the vapor is a non ideal mixture. They are 
given by 



Rn2 " 42^N2 <''N2 - %2”^ 

with a similar expression for P 



O 2 



Pn 2>/RT - f 6 n202<='02>^/'’'^ 

and where 



( 5 n 202 ^%202 " 



(69) 



If 

= 



^N 2 ^N 2 ^ 02^62 



^2 ’'02 a + b(xj^^ 



X 



02 



) + c(xL-xJ^)2. . . 



N2 ^02' 



(70) 



then 
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(71) 



and similarly for In'^Q^. 

The values of 1, ra and n are expressed as 




(72) 




and A = a/RT, B=b/RT, C=c/RT . 

As a first approximation the solution is assumed to be regular, 

such that B=C=0, By neglecting corrections for the vapor phase 

non ideality, A is set equal to 4ln(2P /P. + P,, ), 

x=.5 02 l>'2 

P,, and P are then calculated for the various experimental 

Nj 02 

liquid compositions from Eqs, (68). These are used to calculate 



with similar expressions for dP/dB and dP/dC using m and n re- 
spectively. 

Changes in A, B & C to most nearly reduce the pressure residuals 
to zero are determined by fitting by least squares to 




and P(>alc Eqs. (71) & (67) which in turn yield 



the pressure residual 



R = P 



- P 

expt calc 



(73) 



and the derivatives 



dP/dA 




(74) 



(dP/dA)6A + (dP/d3)6B + (dP/dC) 6c = R 



(75) 



05 



or by solving simultaneously equations of the form 

(5A^(dP/dA)^ + 6 b ^(dP/dA)(dP/dB) + (dP/dA) (dP/dC) 

^R(dP/dA) 

with similar expressions for ^R(dP/dB) and ^R(dP/dC), 



( 76 ) 



Appendix C 



Selected Data 



1. Squilibrium Data 



p 


T 


xl 








si 


S^ 


Atm 




Oracle ^2 


Oracle N 2 


J/mole 


J/niole 


j/mole °K 


j/mole 


1 


77.36 


100 


100 


823.7 


6793.5 


92.84 


165.00 




78 


89.8 


97.3 


722.0 


6790.0 


95.62 


166.23 




80 


64.0 


87.6 


472.0 


6727.0 


98.00 


169.18 




82 


44.3 


77.0 


283.0 


6666 . 0 


98.27 


171.28 




84 


29.6 


63.8 


160.0 


6607.0 


97.88 


173.05 




86 


18.0 


47.6 


70.0 


6547.0 


97.08 


174.00 




88 


8.3 


28.0 


8.0 


6486 . 0 


95.95 


173.88 




90 


1.0 


3.2 


-29.0 


6426 . 0 


94.50 


170.73 




90.18 


0 


0 


-32.0 


64 o 6 . 1 


93.94 


169.63 


5 


94.18 


100 


100 


1811.5 


6672.3 


104.21 


155.84 




96 


79.3 


91.0 


1630.0 


6735.0 


108.25 


159.05 




98 


61.3 


80.4 


1472. 0 


6800.0 


109.49 


161.33 




100 


45.7 


■69.2 


1337.0 


6865 . 0 


109.60 


163.09 




102 


33.0 


56.? 


1245.0 


6928.0 


109.28 


164.32 




104 


22.0 


42.4 


1171.0 


6994 . 0 


108.58 


164.80 




106 


12.3 


26 . 6 


1108.0 


7053.0 


107.45 


164.58 




108 


3.6 


9.3 


1057.0 


7116.0 


105.78 


162.90 




108.99 


0 


0 


1037.9 


7144.7 


104.60 


160.63 


10 


103.95 


100 


100 


2440.2 


6702.3 


110.34 


151.36 




104 


99.3 


99.5 


2433.0 


6704.0 


110.57 


151.53 




106 


79.3 


89.0 


2265.0 


6783.0 


114.27 


■155.12 




108 


63.6 


78.7 


2137.0 


6855.0 


115.38 


157.24 




110 


49.6 


68.0 


2028.0 


6926.0 


115.65 


158.88 




112 


37.3 


56.2 


1939.0 


6994.0 


115.^9 


159.95 




114 


26 o 3 


43.4 


1858.0 


7060.0 


114.89 


160.55 




116 


16.2 


29.1 


1784.0 


7125.0 


113.79 


160.50 




118 


6.9 


13.8 


1721.0 


7186.0 


112.12 


159.30 




119.86 


0 


0 


1681.7 


7245.1 


110.09 


156.50 


20 


115.82 


100 


100 


3321.9 


6515.4 


117.96 


145.53 




116 


98 ,'2 


99.0 


3315.0 


6529.0 


118.48 


146.00 




118 


82.3 


89.1 


3177.0 


6649.0 


121.52 


149.78 




120 


68.8 


79.4 


3072.0 


6752.0 


122.79 


151.98 




122 


56.1 


69.2 


2972.0 


6839.0 


123.09 


153.70 




124 ■ 


44.1 


58.0 


2873.0 


6919.0 


122.98 


154.91 




126 


33.0 


46.0 


2785.0 


6993.0 


122.46 


155.62 




128 


22.8 


33.7 


2707.0 


7058.0 


121.53 


155.77 




130 


13.1 


20.7 


2631.0 


7125.0 


120.13 


155.23 




132 


4.1 


7.3 


2543.0 


7190.0 


118.13 


153.77 




133.03 


0 


0 


2542.0 


7208. 3 


116.61 


151.69 
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2. Pure Component Data 



P 

Atm 

1 



5 



10 



20 



T 


HO 2 


8N2 


S02 o 


2^2 ^ 


°K 


J/mole 


j/moie 


j/mole °K 


J/raole °K 


65 


-1347.1 


108.4 


76.95 


82.77 


70 


-1115.3 


397.6 


80.38 


87.05 


75 


-861.6 


687.0 


83.88 


91.05 


80 


-594.3 


6486.3 


87.32 


166.02 


82 


- 485.3 


6546 . 8 


88.67 


166.77 


8^ 


-375.9 


6607.3 


89.99 


167.50 


86 


-265.2 


6667 . 8 


91.30 


168.21 


88 


-154.0 


6728.0 


92.57 


168.90 


90 


- 42.3 


6788.3 


93.83 


169.56 


95 


6940.0 


6938.4 


171.21 


171.20 


100 


7091.1 


7088. 1 


172.76 


172.73 


80 


- 587.5 


983.1 . 


87.28 


94 . 68 


85 


-313.5 


1270.8 


90.60 


98.18 


90 


-135.6 


1561.9 


93.78 


101.52 


96 


300.9 


6737.0 


97.40 


156.51 


100 


526.4 


6876.3 


99.70 


157.91 


104 


752.9 


7012.4 


101.92 


159.26 


108 


981.0 


• 7145.8 


104.07 


160.52 


110 


7179.5 


7211.9 


160.94 


161.13 


115 


7349.4 


7374.9 


162.45 


162.56 


120 


7515.3 


7535.7 


163.87 


163.93 


100 


533.9 


2178.0 


99.63 


107.79 


104 


759.8 


6704.5 


101.84 


151.38 


106 


873.2 


6790.2 


102.92 


152.20 


108 


987.1 


6872.9 


103.99 


•152.95 


110 


1101.6 


6953.3 


105.04 


153.71 


112 


1216.9 


7031.7 


106.07 


154.41 


114 


1333.1 


7108.5 


107.10 


155.08 


116 


1450.5 


7183.6 


• 108.12 


155.75 


118 


1569.5 


7257.5 


109.14 


156.37 


120 


7250.7 


7330.7 


156.55 


156.99 


105 


830.3 


2508.0 


102.23 


110.59 


110 


1113.5 


2856.0 


104.86 


113.84 


115 


1400.8 


3249.0 


107.41 


117.32 


118 


1576.5 


6663 . 6 


108.92 


146.89 


122 


1816.7 


6894.7 


110.92 


148.72 


126 


2066 .4 


7096.5 


112.94 


150.35 


130 


2329.7 


7280.5 


115.00 


151.80 


135 


7312.6 


7494.8 


152.46 


153.40 


140 


7552.3 


7697.1 


154.21 


15^.89 


145 


7770.4 


7890.4 


155.74 


156.23 



I ' 

HI 




■f 

I 
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3. Deviations 





X 


'C' 


s'*' 


S^ 




^mole 


j/roole 


j/mole °K j/mole °K 




0 


0 


0 0 






10 


36.4 


2.73 


26 




20 


56.1 


4.14 


35 




30 


65.7 


5.09 


38 




40 


68.6 


5.60 


37 




50 


66 . 9 


5.73 


3 ^ 




60 


62.3 


5.60 


29 




70 


55.2 


5.09 


28 . 




80 


^ 3.9 


4.14 


24 




90 


26.8 


2.73 


18 




100 


0 


0 0 




4, Summary of equimolar G 






oT 


gE 


Fig . 33 


Source of 


Source of vapor 




j/mole 


notation 


calculation 


pressure data 


63.14 


10.88 


A 


7 


7 


65 


10.76 


B 


1 


9 


70 


10,28 


It 


tt 


tt 


77.5 


9 . 7 ^ 


It 


It 


tt 


75.16 


9.65 


c 


5 


.8 


77.5 


10.75 


D 


II 


5 


77 


10.0 


E 


2 


2 


83.82 


9.2 


F 


6 


6 


90.0 


9.6 


G 


It 


10 


91.0 


7.3 


tt 


II 


10 


80 


8c45 


H 


This work 


8 


90 


8.82 


tt 


tt 


tt 


100 


9.75 


ti 


tt 


tt 


110 


9.08 


tt 


ft 


tt 


120 


7.64 


tt 


tt 


tt' 



Source notation refers to "References for Excess Functions" 
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